1 Purpose

This document summarizes updated methods and results from the analyses origionally described here. Notably, our analysis now includes Pink (n=70) and Chum (n=44) salmon stocks from across the eastern North Pacific, as well as updated data for Sockeye stocks (n = 52).

2 Methods

The overarching goal of our analyses is to characterize relationships between ocean conditions (SST and abundance of potential competitors) and salmon productivity (R/S) over space and time. Since our last update, we have made minor changes to our approach, as described in the sections below.

2.1 Data

We have compiled spawner-recruitment time series for all five species of Pacific salmon, which are hosted in a github repository, although only Sockeye, Pink, and Chum populations are analyzed here. We analyze a total of 166 populations originating from Washington to the Bering Sea:

  • From our original set of 58 Sockeye populations, we obtained more data for the Southeast Alaska and West Coast (BC) regions, but subsequently chose to remove time series that did not span at least 1984-2015 brood years without significant gaps. This is intended to make sure each time series overlaps with all three time periods in our ‘era’ models, and excludes short series that random walk models may be sensitive to. There are now 52 Sockeye stocks in the analysis - many of the original stocks, and some new ones.

  • We obtained spawner-recruitment time series for Pink and Chum populations, and currently use all time series with reasonable data quality, regardless of whether they meet the criteria applied to Sockeye data. These data are not representative of all regions and time periods, and notably, we have no data after 2009 for any Alaskan stocks. We are working with Andrew Munro (ADF&G) to update these Alaskan timeseries (or some subset of them) through to the present.

  • Compilation of Coho and Chinook time series is on the back burner for now.

2.1.1 Map of ocean entry points

The map below shows the ocean entry locations of all stocks in the analysis presented. The legend to the right of the plot allows you to toggle between species (double-click a legend item to isolate that species). Hovering over a point will show the stock name and ocean region.

Figure 1. Ocean entry locations of Sockeye (n=52), Pink (n=70), and Chum (n=44) salmon stocks included in analyses (total N=166).

2.1.2 Time Series Length

Productivity (log[R/S]) time series of all stocks are shown below to illustrate the length and relative number of time series among regions and species. Vertical dashed lines represent proposed regime shifts in the ‘era’ model (described in sec. 2.2) to illustrate the coverage of data over periods of interest. Where coloured lines break to light gray, data are missing.

2.1.2.1 Sockeye

Figure 2. Productivity (log[R/S]) time series of Sockeye stocks (n=52) with vertical dashed lines indicating proposed ocean regime shifts.

2.1.2.2 Pink

Figure 3. Productivity (log[R/S]) time series of Pink stocks (n=70) with vertical dashed lines indicating proposed ocean regime shifts. Note that Even- and Odd-run pinks are considered separately due to their independent and sometimes vastly different productivity regimes.

2.1.2.3 Chum

Figure 4. Productivity (log[R/S]) time series of Chum stocks (n=44) with vertical dashed lines indicating proposed ocean regime shifts. Note that there are no time series in the Southeast Alaska (SEAK) region; one stock (Middle Skeena) that would have been assigned to SEAK was grouped with West Coast stocks instead.

2.2 Models

Our updated analyses focus on three classes of generalized spawner-recruitment models:

  1. Stationary models (e.g. Connors et al. 2020), which estimate time-invariant relationships. These have not changed since last update other than the addition of data for the two new species.

  2. ‘Era’ models (e.g., Malick 2020), which allow relationships to vary among pre-defined periods that represent hypothesized shifts in North Pacific Ocean processes and relationships.

    • We no longer model 1976/77 as a regime shift as per feedback from our last update, which highlighted that this shift from a cold to warm PDO phase simply reflected a temperature change (which is already accounted for via our direct consideration of ocean temperature) and not a fundamental change in oceanographic processes like the 1988/89 regime shift.

    • Our updated analyses consider a potential regime shift at the onset of the ~2013-2015 marine heatwave (‘the blob’). Currently, we consider brood years >= 2011 to have interacted with the marine heatwave regardless of the species and stock. In future analyses, the timing can be altered to better reflect the diverse life histories of populations in the analysis.

  3. Random walk models (e.g., Malick 2020), which allow relationships to evolve gradually through time.

    • These have not changed since last update other than the addition of new data. However, a version that estimates ocean basin-scale trends in relationships (e.g., as shared trends within ocean region that are modeled hierarchically) is currently in initial stages of development.
  1. At this time we are not pursuing Hidden Markov models, which model the system as having two or more latent states, as our initial analyses indicated they did not add any additional inference beyond those from the era and random walk models. They may be revisited in the future.

Further details on each model class and Bayesian model fitting procedure can be found here.

3 Results

3.1 Overall insights

There is some evidence for spatially and temporally varying responses to warming and competition at sea in all three species analyzed.

In Sockeye, we continue to see strong evidence of spatial variation in relationship between productivity and both SST and competitor abundance, as has been described previously, as well as variation in relationships within regions over time. For southernmost stocks we see a general strengthening of inverse relationships for both covariates over time. For stocks at intermediate latitudes (SEAK and GoA) we see evidence of relationships with SST flipping over time such that during the ‘marine heatwave’ period they are strongly negative, and we see weak and variable relationships in early and late periods and consistently negative ones in the 1989-2010 period. At the northernmost extent (Bering Sea) we see consistently positive relationships with SST and variable relationships with competitor abundance. Due to the short duration of the marine heatwave period that we are able to consider, we should interpret findings for the period cautiously.

In Pink salmon, we see clear evidence of spatial variation in SST relationships, consistent with previous work, but less pronounced patterns in ones with competitor abundance. Like in Sockeye, relationships with SST vary latitudinally and for southernmost stocks become more negative over time while for SEAK stocks there is evidence of a flip in relationships from positive to negative by the ‘marine heatwave’ period. For the competitor index we see positive relationships in the early period which presumably reflects shared positive responses to ocean conditions and then highly uncertain and weak relationships thereafter which is consistent with the hypothesis that Pink salmon are the least affected by intra-specfic competition at sea (Ruggerone & Nielsen, 2005). However, in addition to the short temporal window of the most recent period, many of the time series for mid and northern latitude stocks have not been updated yet.

In Chum, there is a consistently negative relationship with SST in the southernmost (WC) region that decreases over time, while more northern stocks show weakly positive to neutral SST relationships. Chum have relatively weak relationships with competitor abundance compared to Sockeye, but northern (BS) stocks shift from a positive relationship in the first era to weakly negative in the second era in both species. Again, this likely represents a response to favourable ocean conditions experienced by all species in the early period, and negative responses to inter-specific competition in later periods as pink salmon abundance increased. Inference about Chum in the marine heatwave period will also require updating time series for Alaskan stocks.

3.2 Stationary models

3.2.1 Sockeye

Figure 5. Posterior probability distributions of the predicted effect of SST (top), competitor (pink salmon) abundance (middle), and the combined effect (bottom) on Sockeye productivity (R/S). Faint lines show stock-specific effects while bold lines show regional effects from hierarchical model. X-axis values represent the percent change in productivity per standard deviation unit increase in the covariate.

3.2.2 Pink

Figure 6. Posterior probability distributions of the predicted effect of SST (top), intra-specific competitor abundance (middle), and the combined effect (bottom) on Pink salmon productivity (R/S). Faint lines show stock-specific effects while bold lines show regional effects from hierarchical model. X-axis values represent the percent change in productivity per standard deviation unit increase in the covariate.

3.2.3 Chum

Figure 7. Posterior probability distributions of the predicted effect of SST (top), competitor (pink salmon) abundance (middle), and the combined effect (bottom) on Chum salmon productivity (R/S). Faint lines show stock-specific effects while bold lines show regional effects from hierarchical model. X-axis values represent the percent change in productivity per standard deviation unit increase in the covariate.

3.3 Era models

3.3.1 Sockeye

Figure 8. Posterior probability distributions of the predicted effect of SST and competitors on Sockeye productivity over three pre-defined time periods/eras (earliest in top panel). Regional mean effects are shown by bold lines and individual stocks’ distributions by light lines.

3.3.2 Pink

Figure 9. Posterior probability distributions of the predicted effect of SST and competitors on Pink salmon productivity over three pre-defined time periods/eras (earliest in top panel). Regional mean effects are shown by bold lines and individual stocks’ distributions by light lines. Note that data from Alaska for the most recent time period are needed (Figure 3).

3.3.3 Chum

Figure 10. Posterior probability distributions of the predicted effect of SST and competitors on Chum salmon productivity over three pre-defined time periods/eras (earliest in top panel). Regional mean effects are shown by bold lines and individual stocks’ distributions by light lines. Note that data from Alaska for the most recent time period are needed, and there are no time series in the Southeast Alaska region.

3.4 Random Walk models

3.4.1 Sockeye

Figure 11. Time-varying posterior mean estimates of SST and Competitor covariate effects on Sockeye productivity, modelled as a random walk. Individual stock estimates are in faint lines, while regional means are represented by bolded lines. Region-wide means are post-hoc calculations, rather than resulting from a hierarchical model structure as in the stationary and era models.

3.4.2 Pink

Figure 12. Time-varying posterior mean estimates of SST and Competitor covariate effects on Pink salmon productivity, modelled as a random walk. Individual stock estimates are in faint lines, while regional means are represented by bold lines. Solid bold lines represent odd-year Pink stocks, while dashed lines are even-year stocks. Region-wide means are post-hoc calculations, rather than resulting from a hierarchical model structure as in the stationary and era models.

3.4.3 Chum

Figure 13. Time-varying posterior mean estimates of SST and Competitor covariate effects on Chum salmon productivity, modelled as a random walk. Individual stock estimates are in faint lines, while regional means are represented by bolded lines. Region-wide means are post-hoc calculations, rather than resulting from a hierarchical model structure as in the stationary and era models.

4 Next steps

Priorities for future analysis include:

  1. Updating data to present (as recent as possible) for all species and regions, particularly Alaskan pink and chum.

  2. Sensitivity analyses of SST and competitor index covariates.

    • SST: The resolution of the SST data is currently 2 x 2 degrees (ERSST). Cells within 400 km of ocean entry points are averaged over the three months after migration to sea (April-September, depending on ocean region according to Mueter et al. 2002). Given that these assumptions about the spatial and temporal window over which SST is averaged are quite dated, sensitivity to these assumptions will be tested with alternative spatial and temporal windows. However, we will continue to focus on early marine life to remain consistent with the overall hypothesis that juvenile salmon are impacted by temperature in the first year of marine life, and fine-scale analyses of local SST conditions are outside the scope of this project. That said there may be high-resolution SST data available in the near future to further explore this in subsequent paper.

    • Competitor Index: Our base cases analysis focuses on the abundance of pink salmon across the North Pacific in the second and/or third years of marine life as an index of potential direct or indirect competition for food. This approach is consistent with research that has suggested some salmon species primarily exhibit responses to competition with pink salmon during their second and/or third growing seasons at sea (Connors et al. 2020; Ruggerone et al.2023). However, we also plan to fit models to a small subset of alternative competitor indices as part of sensitivity analysis for this paper (e.g., total pink, chum and sockeye abundance; just North American abundances, biomass instead of abundance, etc.) and have begun sketching out a more comprehensive analytical plan for a separate paper that explores evidence for a wider range of hypotheses about intra- and inter-specific interactions at a variety of spatial and temporal scales.

  3. Revising hierachical grouping of Pink salmon stocks. The random walk models for Pink salmon show the trends for even-year and odd-year stocks separately for ease of visualization, but the difference in the effects between lineages in some cases begs the question of whether even and odd stocks within a region should be grouped together or separately in hierarchical models. We will investigate the possibility of separating them, which is consistent with the hypothesis that the genetically distinct odd- and even-year lineages of Pink salmon may respond differently to ocean warming (Irvine et al. 2014).

  4. Random Walk model improvements. Work on versions of random walk models with (1) estimated region-wide mean effects, and (2) reduced sensitivity to time series of different lengths, is ongoing.